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We present results for the large-pT inclusive tt" production in p-p and A-A collisions at RHIC and 
LHC energies. We include the full next-to-leading order radiative corrections, 0{al), and nuclear 
effects such as parton energy loss and nuclear shadowing. We find the next-to-leading order correc- 
tions and the parton energy loss effect to be large and pr dependent, while the nuclear shadowing 
effects are small (< 10%). We calculate the ratio of prompt photons to neutral pions produced in 
heavy ion collisions and show that at RHIC energies this ratio increases with pr approaching one 
at pt ~ 8 GeV, due to the large suppression of 7r° production. We show that at the LHC, this 
ratio has a steep pr dependence and approaches 25% effect at pr ~ 40 GeV. We discuss theoretical 
uncertainties inherent in our calculation, such as choice of the renormalization, factorization and 
fragmentation scales and the K-factors which signify the size of higher-order corrections. 



I. INTRODUCTION 

In high energy heavy ion collisions hard scatterings of 
partons occur in the early stages of the reaction, well 
before the possible formation of a quark gluon plasma, 
resulting in production of large transverse momentum 
particles. Fast partons produced in these hard collisions 
propagate through the hot and dense medium created in 
the heavy ion collision and lose their energy. Therefore, 
a possible signature for parton energy loss in a hot, dense 
medium is the suppression of pion production in heavy 
ion collisions relative to hadron-hadron collisions. The 
Relativistic Heavy-Ion Collider (RHIC) at BNL with Au- 
Au colhsions at = 130 GeV and at = 200 GeV, 
and the Large Hadron Collider (LHC) at CERN which 
will coUide Pb-Pb at -^s = 5.5TeV provide the best op- 
portunity to study the properties of the hot and dense 
matter and the possible formation of a new phase, the 
quark-gluon plasma. Recent measurements of inclusive 
7r° production at RHIC energy of ^/s = 130 GeV [1] and 
at ^/s = 200 GeV [2] show a large suppression in the pr 
region (lGeV< pt < lOGeV) which has created much 
excitement in the field. The question of the origin of this 
suppression has inspired many interpretations. Clearly, 
the low Pt region has large nonperturbative contributions 
and perturbative calculation cannot be trusted. Here we 
concentrate on production of tt" in the large-pr region 
where perturbative QCD calculations are expected to be 
more reliable. 

In addition to being of interest for studying nuclear ef- 
fects, such as parton energy loss and nuclear shadowing, 
large-pr mesons form a significant background for the 
prompt photons. In principle, large-pT pions could form 
"fake prompt photons" when one photon from pion decay 
escapes detection. Theoretical predictions for the ratio of 
prompt protons to pions at RHIC and LHC energies are 
crucial for studying possible quark-gluon plasma forma- 
tion via photons. In addition, this ratio may reduce some 



of the theoretical uncertainties, such as choice of factor- 
ization, renormalization and fragmentation scale, or the 
choice of gluon fragmentation function. 

In perturbative QCD, the inclusive cross section for 
pion production in a hadronic collision is given by: 



U-^Ptv 



dxa / dx 



dz^F,{Xa,Q^) 



Fj{xb,Q^)D,/Jz,Q})E^ 



d^a,. 



d^PTT 



(1) 



where Fi{x,Q^) is the i-th parton distribution in a nu- 
cleon, Xa and Xb are the fractional momenta of incoming 
partons, Dc/Tr{z,Q'j) is the pion fragmentation function, 
z is the fraction of parton energy carried by the pion and 

are parton-parton cross sections which include 



leading-order, 0(af), subprocesses such as: 
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and the next-to-leading order, 0{a^), subprocesses such 
as: 
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The running coupling constant UsifJ.'^), calculated to 
next-to-leading order, is given by 
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where ^ is the renormahzation scale, Np is the number 
of flavors and A is the Aqcd scale. 

The parton distribution functions Fi{x,Q^) are mea- 
sured in Deep Inelastic Scattering experiments such 
as those at HERA [3] while fragmentation functions, 
Dc/-n{z, Q/), that describe the transition of the partons 
into the final-state pions are extracted from e^e~ an- 
nihilation data from PETRA, PEP and LEP [4]. The 
gluon fragmentation function, which gives the dominant 
contribution to tt*^ production at the LHC, is not well 
determined by e+e~ data, since it appears only at NLO. 
Nevertheless, it is possible to get some constraint on the 
gluon fragmentation function from measurements of large 
Pt pion production in hadronic collisions at high energies. 
Such a study has been done using UAl data in the range 
5GeV< Pt < 20GeV [5]. The gluon fragmentation func- 
tion of Binnewies, Kniehl and Kramer (BKK) [4] is found 
to be consistent with the UAl data. The fragmentation 
scale, Qf, when taken to be too small, i.e. Qf = Pt/^i 
probes the region currently not tested by the data. In ad- 
dition, theoretical improvement is needed in resumming 
large ln{l — z) terms present in the higher-order correc- 
tions. 

Divergencies present in the calculation of parton cross 
section, daij, require careful treatment of coUinear and 
"soft" singularities in the matrix elements. Divergences 
cancel between real and virtual graphs when the physical 
7r° cross section is calculated. Infrared or "soft" diver- 
gences {k — > 0) cancel between virtual amplitudes and 
real amplitudes. Collinear divergences {9 ^ i) between 
initial quark and radiated gluon, for example) are reab- 
sorbed in structure function and fragmentation function 
renormahzation ("mass" singularities). 

Detailed study of large-pr inclusive tt" production in 
hadronic collisions show very good agreement between 
theory and experiments, apart from an overall normal- 
ization [6]. 



II. large-Pt inclusive tt" production in 

HEAVY ION COLLISIONS 

To calculate the inclusive cross section for pion pro- 
duction in heavy ion collisions, we will use (1) with the 
distribution and fragmentation functions appropriately 
modified to include nuclear effects such as shadowing and 
energy loss. 

It is a well known experimental fact that the distri- 
bution of quarks and gluons inside nuclei is modified 
compared to that in a free nucleon. This modification 
is known as nuclear shadowing (for a review of nuclear 
shadowing, we refer the reader to [7]). The parton dis- 
tribution in a nucleus Fa/A{x, Q^, bt), can be written as 

F,/a{x, Q\bt) = TA{ht) S,,a{x, g2) F,/^{x, g2) 



where TA{bt) is the nuclear thickness function, 
^a/Ni^jQ'^) is the parton distribution function in a nu- 
cleon and 5*^/^(2;, Q^) is the parton shadowing function, 
Sa/A{x,Q^) = Fa/A{x,Q^)/AF,/N{x,Q^). It should be 
emphasized that the dependence of nuclear shad- 
owing is poorly known, especially in the small Xbj re- 
gion. Also, gluon shadowing is measured only indirectly 
through scaling violation of F2 structure function which 
leads to large uncertainties in the gluon distribution func- 
tion in nuclei. In this work, we use the nuclear shadow- 
ing parametrization due to Eskola, Kolhinen and Sal- 
gado (EKS98), which is dependent and distinguishes 
between quarks and gluons [9] and was shown to be in 
very good agreement with the NMC data on depen- 
dence of F/"/F2^ [10], while some other parametriza- 
tions which have large gluon suppression due to nuclear 
shadowing are ruled out [11]. It is also shown that mod- 
ifications to the DGLAP evolution due to gluon fusion 
are small for the kinematic region of relevance to RHIC. 
In Fig. (1) we show EKS98 nuclear shadowing function 
at — 2.25GeV^ and, for comparison, we also show 
nuclear shadowing parametrization due to Benesh, Qiu 
and Vary (BQV) [8] which is independent and treats 
shadowing of quarks and gluons on the same footing. 
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FIG. 1. The nuclear shadowing ratio as parametrized by 
BQV and EKS98. 

In Fig. (2) wc show the dependence of the shad- 
owing function for gluons in EKS parametrization [12]. 
Clearly, there is a strong dependence in the EKS 
parametrization of nuclear shadowing, especially for glu- 
ons, as well as large anti-shadowing effect in the large Xbj 
region. 
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Another nuclear effect which we include is the medium 
induced energy loss. It is expected that fast partons prop- 
agating through the hot and dense medium created after 
a high energy heavy ion collision will scatter from the 
partons in the medium, lose part of their energy and 
then fragment into hadrons with a reduced energy. As a 
result, the spectrum of the final state hadrons observed 
in a heavy ion collision is expected to be suppressed com- 
pared to hadronic collisions. This has been extensively 
studied and we refer the reader to [13] for some recent 
work. The Landau-Pomeranchuk-Migdal (LPM) effect 
occurs when an ultrarelativistic particle emits low energy 
bremsstrahlung photons (gluons) as the particle passes 
through dense matter; fewer photons (gluons) are emit- 
ted than predicted by bremsstrahlung theory for isolated 
atoms. 

While a dynamical study of the parton propagation 
in a hot and dense medium created in a realistic heavy 
ion collision and the modification of the hadronization 
is more desirable, we will use a phenomenological model 
[14] here to demonstrate how sensitive pion production 
is to the average energy loss suffered by a parton in a hot 
and dense medium. We restrict ourselves to the central 
rapidity region so that a parton will only propagate in the 
transverse direction in a cylindrical system. Given the 
inelastic scattering mean- free-path, Aa, the probability 
for a parton to scatter n times within a distance AL 
before it escapes the system is assumed to be given by 
the Poisson distribution 
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In the Huang-Sarcevic-Wang model of energy loss, the 
hadronic fragmentation function Dc/Tri^, Q^) is modified 
in order to include multiple scattering of a parton in a 
nuclear medium. If the average energy loss per scattering 
suffered by the parton a is ea, the nuclear fragmentation 
functions can be modeled as [14], 
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where z° — z/[\ — nea/ET)^ z'^ = zETj^a, N is the 
maximum number of collisions for which z°: < 1 and , 

n — C/TT 

is the hadronic fragmentation function which gives the 
probability that quark or a gluon would fragment into a 
pion. The first term corresponds to the fragmentation of 
the leading partons with reduced energy Ex — nta and the 
second term comes from the emitted gluons each having 
energy ta on the average. The average number of scat- 
terings within a distance Ai is (ria) = AL/Xa- We take 
Aa = ifm and AL = Ra- 

In this work, we use the NLO code due to Aurenche 
et al. for calculating tt*^ and prompt photon produc- 
tion in p-p collisions [6] with MRS99 parametrization 
of nucleon structure functions [15] and BKK pion frag- 
mentation functions [4] . We use EKS98 shadowing func- 
tions to include nuclear shadowing and we modify the 
BKK fragmentation functions according to the model 
of Huang-Sarcevic-Wang [14] in order to take into ac- 
count medium induced parton [6] energy loss effects. We 
calculate the invariant cross section for tt" and prompt 
photon production in heavy-ion collision normalized to 
the number of binary nucleon- nucleon collisions, Ncoiu 
where Ncoii can theoretically be determined from nu- 
clear overlapping function, i.e. Ncoii — TAA(b)(J^^i and 

TAAib) = J (fbiTAil h DTaH b- bi I). The number of 
N-N collisions depends on the centrality that experiment 
triggers on. Here we take NcoU = 975, which is obtained 
by PHENIX for their central collisions [2]. In order to 
investigate the sensitivity of our results to the choice of 
energy loss parameters in the model of [14], we consider 
various forms of the average parton energy loss parame- 
ter, e; constant energy loss as well as two different energy 
dependent forms of energy loss. We find that fractional 
energy loss, ^ ~ 6.1% gives the best description of the 
ttO data at RHIC for = 2mGeV (see Fig. 8). With 
this energy loss, we make predictions for prompt photon 
invariant cross sections. We set all the scales, factor- 
ization, renormalization and fragmentation scale to be 
equal, Q = Qp = fi = npx where we take k = 1 and 
K = 2. We discuss sensitivity of our results to the choice 
of scale* We assume that nucleus consists of A protons, 
i.e. we do not take into account isospin effects, which are 
expected to be small for large A. 

In Fig. (3) we show the invariant cross section -^r^ 
for inclusive tt*^ production in proton-proton and nucleus- 
nucleus {A = 200) collisions at = 130 GeV. There 



*We do not consider k = 0.5 because for low pt the scales 
would be in the region where currently there is no data and 
where perturbative QCD in not reliable. 
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is a clear suppression of the nuclear cross section com- 
pared to the hadronic one. While constant energy loss 
of e = O.SGeV and LPM type energy loss e ~ \/E show 
similar shape of pt spectrum, fractional energy loss (BH) 
leads to a much steeper spectrum. We choose values of 
e and Elpm such that the suppresion at px = 4GeV 
is in agreement with the data. However, the data at 
= 130Gey is not precise enough to distinguish be- 
tween the different forms of energy loss. Nuclear shad- 
owing effect, which was incorporated using EKS98 shad- 
owing, is found to be very small (few percent) at RHIC 
energies. 
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FIG. 3. Inclusive 7r° distribution at ^ = 130 GeV. 

From Fig. (3) we see that the 7r° production cross 
section in heavy ion collisions is suppressed compared to 
hadronic collisions even at very large pr {pt > 10 GeV). 
This is due to the fact that high pT partons lose their 
energy as they propagate through the dense medium and 
contribute to the lower pr pions. 

The next-to-leading order contributions, O(a^), which 
are included in our calculations, are large and pT depen- 
dent. For example, the nuclear K- factor defined as a ratio 
of full next-to-leading order calculation to the leading or- 
der is varying between 1.6 and 2.2 for px between 4GeV 
and 20GeV, independent of the form of parton energy 
loss used. This is shown in Fig. (4). One should keep 
in mind that next-to-leading order calculation in nuclear 
collisions is not as complete as for the pp case due to 
the fact that EKS nuclear shadowing function is based 
on the leading-order DGLAP evolution and the parton 
energy loss has been calculated only in the leading-order 
approximation. Uncertainties in the modified fragmen- 
tation function are most likely much larger than the un- 
certainty due to the next-to-leading order corrections to 
the nuclear shadowing function. 
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FIG. 4. K-factor, defined as the ratio of full NLO and LO 
tt" inclusive cross section at ^/s = 130 GeV. 

In Fig. (5) we show the ratio of inclusive tt" cross sec- 
tion in nucleus-nucleus collisions to the one in p-p colli- 
sions. Since the nuclear shadowing effects are small, most 
of the observed suppression is due to energy loss effects. 
The difference between different forms of energy loss is 
striking, specially in the low to intermediate px region 
where they have different slopes. This difference is even 
more pronounced for the case of fractional energy loss 
where this ratio actually decreases. However, we point 
out that we do not expect our rather crude model of en- 
ergy loss through modified fragmentation functions to be 
correct at very high pT ■ This is due to the fact that with 
a constant fractional energy loss as in the BH case, the 
ratio of AA and pp cross sections will never become one as 
one goes to higher momenta unlike the LPM case where 
fractional energy loss goes away as 1/Ve. Although we 
cannot determine exact value of pT at which our model 
will break down from first principles, it is clear that it 
will happen at higher values of pt for higher energies. 
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FIG. 5. Ratio of inclusive tt" distributions in heavy ion and 
p-p collisions at -y/s — 130 GeV. 
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We calculate prompt photon production in heavy ion 
collisions using Aurenche et al. program for prompt pho- 
ton production in p-p collisions that includes all next-to- 
leading order contributions [6]. In prompt photon pro- 
duction there are two types of subprocesses that con- 
tribute, direct processes and bremsstrahlung processes, 
the later one being convoluted with the fragmentation 
functions, D^/^{z, Q^), that describe the transition of the 
partons into the final-state 7. The fragmentation func- 
tions without medium effects are extracted from e+e" 
data [4]. We use MRS99 parametrization of nucleon 
structure functions [15], EKS98 nuclear shadowing func- 
tions and we modify fragmentation functions in the same 
way as for pions. 

In Fig. (6) we show the inclusive cross section for 
prompt photon production in heavy ion collisions at 
= 130 GeV, normalized to proton-proton collisions, 
together with the same ratio for inclusive tt*^ produc- 
tion. We note that for the same energy loss param- 
eter e — 6.1%E , production is suppressed more 
than prompt photons. This is due to the fact that at 
Pt = 3GeV, about 75% of the photon contributions 
comes from "direct" processes which are not affected by 
parton energy loss, while in case of pions, all subprocesses 
are convoluted with the modified fragmentation function. 
At higher values of pT, direct photons become even more 
dominant at RHIC and their energy loss becomes negli- 
gible (< 10% at Pt = lOGeV^). Nuclear shadowing effect 
is small both in photon and in pion production. 
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FIG. 6. Prompt photon and inclusive vr" cross sections in 
heavy ion collisions normalized to p-p at ^/s = 130 GeV. 



In Fig. (8) , we show the ratio of invariant cross sections 
for AA and pp collisions at ^/s = 200GeV for different 
forms of energy loss. The pt spectrum behaves very dif- 
ferently in the case of fractional energy loss as compared 
to the constant and LPM type energy losses. Compari- 
son to the preliminary PHENIX data [2] seems to favor 
the fractional energy loss scenario. Furthermore, varying 
the scales from pT to 2pT, in case of LPM energy loss, 
gives large uncertainty, about 30%, while in case of BH, 
this uncertainty is much smaller. This can be seen in Fig 
(8). It is also interesting to note that change of scales 
does not affect the shape of the ratio. 




FIG. 7. Inclusive 71° spectrum at ^ = 200 GeV. 

Prompt photon and pion invariant cross sections nor- 
malized to pp collisions are shown together in Fig. 
(9). Clearly, prompt photons are affected much less 
than pions due to direct photons dominating over 
bremsstrahlung photons at high pT- This is contrary 
to expectations from approaches where the suppresion 
of high Pt spectra is due to high gluon density in the 
initial state (color glass condensate). Therefore, measur- 
ing prompt photons would enable one to tell whether the 
observed suppresion is due to initial state or final state 
effects. 



We show our results for inclusive tt^ production at 
RHIC at Vs = 200 in Fig. (7) for different forms of en- 
ergy loss. Again, the difference between different energy 
loss scenarios is significant especially at the highest pT 
where tt*' production differs by an order of magnitude. 
We find that nuclear shadowing effects are small (less 
than 10%) and the suppression is mostly due to energy 
loss effects. 
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. 8. Ratio of inclusive n" cross sections in heavy ion 
p collisions at = 200 GeV, data is from PHENIX 
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FIG. 9. Prompt photon and inclusive vr" cross sections in 
heavy ion collision normalized to p-p at ^/s — 200 GeV. 

Measurements of inclusive pion production at RHIC 
energies for > 4 GeV could therefore provide valu- 
able information about the medium induced parton en- 
ergy loss since nuclear shadowing effects are very small 
(few %) and most of the observed suppression of hadronic 
spectra in heavy ion collisions is due to energy loss. 

In Fig. (10) we show theoretical uncertainty in predict- 
ing inclusive tt" cross section by varying scales from Q = 
PT to Q = 2pT, with the assumption that Q = Qf = 1^- 
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FIG. 10. Theoretical uncertainty in inclusive tt*^ cross sec- 
tion due to choice of scale for Q — pr and Q — 2pT- We 
assume that Q = Qf ~ jJ-- 

We note that uncertainty is about 40% and that, most 
importantly, the shape of the px spectrum is not very 
sensitive to the choice of scale (the uncertainty due to 
scale choices of Q = 2pT and Q = pt/2 are about 100% 
even though ratio of cross sections is still very insensitive 
to the scale choice). Similar conclusion was reached in 
case of tt" production in hadronic collisions [6]. Com- 
parison of PHENIX data for inclusive 7r° production in 
hadronic collisions shows very good agreement with our 
NLO predictions obtained with all scales set to be equal 
to Pt [2] . The ratio of tt*' production in heavy-ion colli- 
sions to the one in proton-proton collision is less sensitive 
to the change of scales, and in case of BH energy loss, the 
uncertainty is only a few percent. 

In Fig. (11) we show the inclusive 7r° cross sections in 
p-p and heavy ion collisions at LHC energy of ^/s = 5.5 
TeV. In case of heavy ion collisions, we show results again 
for three different forms of energy loss parameter e and 
include nuclear shadowing using the EKS98 parametriza- 
tion. In Fig. (12), we present the ratio of 7r° inclusive 
cross section in heavy ion collisions to p-p. Again, BH 
energy loss leads to a very different behavior of this ratio 
as compared to constant and LPM energy loss. While 
this ratio has steeply increasing pT dependence in the 
case of constant and LPM energy loss, it is almost flat in 
a large pT region in the case of BH energy loss. Again, 
nuclear shadowing effects are small due to the large pt's 
involved. 
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FIG. 11. Inclusive tt" cross sections at LHC energy, 
= 5.5 TeV. 
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FIG. 13. Prompt photon and inclusive n'^ cross sections in 
heavy ion colhsions normalized to p-p at -^s = 5.5 TeV. 
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Finally, in Figs. (14, 15, 16) we show the ratio of 
prompt photons to tt'^ as a function of transverse mo- 
mentum for energies =130 GeV, y/s — 200 GeV and 
^/s = 5.5 TeV, including nuclear shadowing and LPM as 
well as BH forms of energy loss. We also show this ratio 
for proton-proton collisions at the same energy. We find 
that in heavy- ion collisions, the px dependence of the 
7/7r° ratio depends strongly on the form of the energy 
loss. For example, in case of LPM energy loss, this ratio 
increases with px slowly, similar to the case in proton- 
proton collisions. However, in case of BH energy loss, i.e. 
for e = 6.1%E GeV, we find that because of the large 
7r° suppression relative to prompt photons at RHIC en- 
ergies, this ratio increases rapidly with px approaching 1 
at ~ 8 GeV. Similarly, at the LHC this ratio increases 
more than order of magnitude, namely from 0.01 at low 
Pt to about 0.25 at large px- This is a much steeper 
increase than in case of proton-proton collisions, which 
increases only by a factor of 8 in the same px range. 



In Fig. (13) we show cross sections for prompt photon 
and inclusive tt*^ production normalized to p-p using the 
BH form of energy loss. Here we note that suppression of 
prompt photons is comparable to the 7r° case because at 
LHC energies, prompt photon production is dominated 
by bremsstrahlung processes (about 60% contribution at 
Px ^ AGeV) which are modified due to the energy loss 
in a similar way to the tt^ case. 
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FIG. 14. Ratio of prompt photon to tt'^ cross sections at 
= 130 GeV. 
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FIG. 15. Ratio of prompt photon to tt'' cross sections at 
^ = 200 GeV. 
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FIG. 16. Ratio of prompt photon to tt'' cross sections at 
^ = 5.5 TeV. 

III. SUMMARY AND CONCLUSIONS 

We have calculated inclusive pion and prompt photon 
production cross sections in heavy ion and hadronic colli- 
sions at RHIC and LHC energies. We have incorrporated 
next-to-leading order contributions, initial state parton 
distribution functions in a nucleus and medium induced 
parton energy loss by modifying the final state pion frag- 
mentation function. 

We have shown that the nuclear K-factor, the ratio of 
NLO to LO pion cross sections, is large and px depen- 
dent. This PT dependence is specially strong at lower en- 
ergies, i.e. RHIC at ^/s = 130 GeV. We find that nuclear 
shadowing effects are quite small, specially at RHIC at all 
Pt's considered. The EKS98 parametrization of nuclear 
shadowing has a strong p^ dependence which affects the 



cross sections through dependence of the factorization 
scale on pt- If one could measure the 7r° spectra at very 
large pT at RHIC to high accuracy (this would require 
an accuracy of better than 5% at pt ^ 20 GeV which 
is quite unlikely), one could in principle rule out some 
models of nuclear shadowing which are independent. 

We have investigated dependence of our results on 
three different forms of energy loss. We considered con- 
stant energy loss as well as LPM and BH type of energy 
loss. The medium induced energy loss effects in pion 
production are large even at yfs — 130 GeV. However, 
since our perturbative calculation is not applicable at 
low Pt, and the RHIC data at ^/s = 12>QGeV is available 
for PT < 4GeV, we can not distinguish between differ- 
ent functional forms of energy loss from RHIC data at 
y/s = 1306*6^. Therefore, we use the existing data at 
RHIC at y/s = 200 GeV to determine the functional form 
and magnitude of the average energy loss per collision, e. 
We find that BH energy loss describes the data the best. 
We then predict the prompt photon production cross sec- 
tion in Au-Au collisions normalized to the p-p case at 
high Pt (4 GeV< pt < 20 GeV). We also predict the in- 
clusive tt'^ production cross section in heavy ion collisions 
normalized to p-p at the LHC energy, ^/s = 5.5, assuming 
BH energy loss. We demonstrate theoretical uncertain- 
ties due to scale dependence of perturbative calculations 
of this ratio are small (< 15% between Q — pt/2 and 
Q = 2pt). 

Finally, we show results for the ratio of prompt photon 
to tt'^ cross sections at RHIC and LHC, of relevance to 
separating different sources of photon production. We 
show that contribution of direct photons is large at RHIC 
(about 75% at pt = ^GeV) but decreases as the energy 
is increasing, becoming 40% at pt = ^GeV (with 60% 
coming from bremsstrahlung processes) at LHC energies. 
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